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ABSTRACT 

We use a simple model of the dynamics of a narrow-eccentric ring, to put some 
constraints on some of the observable properties of the real systems. In this work we 
concentrate on the case of the 'Titan ringlet of Saturn'. 

Our approach is fluid-like, since our description is based on normal modes of oscil- 
lation rather than in individual particle orbits. Thus, the rigid precession of the ring is 
described as a global m = 1 mode, which originates from a standing wave superposed 
on an axisymmetric background. An integral balance condition for the maintenance of 
the m = 1 standing-wave can be set up, in which the differential precession induced by 
the oblateness of the central planet must cancel the contributions of self-gravity, the res- 
onant satellite forcing and collisional effects. We expect that in nearly-circular narrow 
rings dominated by self-gravity, the eccentricity varies linearly across the ring. Thus, 
we take a first order expansion and we derive two integral relationships from the rigid- 
precession condition. These relate the surface density of the ring with the eccentricity 
at the center, the eccentricity gradient and the location of the secular resonance. 

These relationships are applied to the Titan ringlet of Saturn, which has a secular 
resonance with the satellite Titan in which the ring precession period is close to Titan's 
orbital period. In this case, we estimate the mean surface density and the location of 
the secular resonance. 

Subject headings: Planetary rings. Celestial Mechanics 



1. Introduction 

The dynamical mechanism that maintains the apse alignment of the observed narrow-eccentric 
planetary rings is basically governed by self-gravity (Goldreich and Tremaine 1979), which would 
provide the appropriate contribution to counter-act the differential precession induced by the oblate- 
ness of the central planet. However, predictions of the total mass of the ring produced by this model 
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axe, in general, not in good agreement with the inferred mass of observed eccentric rings (Tyler et 
al. 1986, Graps et al. 1995, Goldreich and Porco 1987). This led to the consideration of other 
factors that might play an important role in the dynamics. In particular, at their narrowest point, 
the ring particles are 'close-packed'. In such a situation particle interaction or pressure effects may 
affect the precession of particle orbits. A simple model where the pinch locks the differential preces- 
sion, was introduced by Dermott and Murray (1980). A more global picture, including the effect of 
stresses due to particle interactions and neighboring satellite perturbations, which offered a better 
agreement with the observations, has also been produced by Bordcries et al. (1983). Their dynam- 
ical model is described in terms of mutually interacting streamlines and the satellite interactions 
(see Goldreich & Tremaine 1981) are computed using a resonance-continuum approximation. The 
standard self-gravity model was later revisited by Chiang and Goldreich (2000), who considered 
the effects of collisions near the edges, proposing that a sharp increase of an order of magnitude in 
the surface density should be observed within the last few hundred meters of the ring edges. More 
recently, employing a pressure term that describes close-packing, Mosqueira and Estrada (2002) 
obtained surface-density solutions that agree well with the currently available mass estimates. 

The eccentric precessing-pattern of the ring can be described as being generated by a normal 
mode of oscillation of wave-number m = 1, which can be viewed as a standing wave. The condi- 
tions for the maintenance of steady global m = 1 modes have been considered by Papaloizou &: 
Melita (2004). To describe the ring perturbations and the m = 1 mode we used the Lagrangian- 
displacement of the particle orbits from their unperturbed circular ones (see for example Shu et 
al. 1985). This model includes the dissipation due to inter-particle collisions, which would lead to 
damping of the mode. However, this global m = 1 mode can also be perturbed by neighboring- 
shepherd satellites, which can inject energy and angular momentum through resonances. In this 
way, losses due to particle collisions can be balanced. Two conditions for the maintenance of the 
rings can be derived. The first one is a condition for the steady maintenance of the amplitude 
or eccentricity associated with the m = 1 mode, which requires the external satellite torque to 
balance the dissipative effects due to collisions (see Papaloizou & Melita 2004). The second one is 
the condition of uniform precession of the ring, which, in the absence of satellite resonances with 
the m = 1 mode, only involves self-gravity and the effect of collisions. These conditions can be 
regarded as continuum forms of the discrete relationships that can be obtained from the 'many 
streamlines' model (Goldreich and Tremaine 1979, Longaretti and Rapapport 1995). 

In this work we produce an extension of Papaloizou and Melita (2004) to the case where 
the pattern frequency of the narrow-eccentric ringlet is in a secular resonance with an external 
satellite. In this case, there is a contribution to the condition of uniform precession, arising from 
the resonant secular perturbation. Resonantly-forced rings are of particular interest, because there 
is a real system that the model can be applied to. It is known that the precession frequency of 
the eccentric Saturnian ringlet at 1.29 Rs (popularly known as the 'Titan ringlet') is in a 1 : 
resonance with the orbital frequency of the Saturnian satellite Titan (Porco et al. 1984). 

From the rigid precession condition two useful relations are derived. If the eccentricity gradient 
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is approximately constant across the ring, which is expected to be the case for a narrow ringlet 

dominated by self-gravity, then, these relationships can constrain the mean surface density, the 
central eccentricity, the eccentricity gradient and the location of the secular resonance. We estimate 
the mean surface density as a function of the location of the resonance and the form of the ring in 
the case where physical collisions are neglected. 

This article is organized as follows. In section 2 we set up the equations for the Lagrangian 
variations starting from the equations of motion in a 2D flat disk approximation. In section 3 we 
give adequate approximations for the evolution of the m = 1 mode when the precession time is 
much longer than the orbital period. In section 4 we derive the condition of rigid precession which 
incorporates secular satellite resonances and we also compute the contribution from the self gravity 
of the ring. Two integral relationships are obtained form the rigid-precession condition in section 5. 
The relation between the eccentricity gradient, the value of the central eccentricity, the location of 
the resonance and the surface density is given in section 6, where we also discuss how to obtain 
the eccentricity gradient from the other parameters in the linear and the non-linear regimes. In 
section 7 we apply these results to the Titan ringlet of Saturn to produce estimates of its mean 
surface-density. Finally, a discussion of the results is given in section 8. 



2. Equations of Motion and Lagrangian Displacement 

We start from the basic equations of motion for a particle in Lagrangian form in 2D: 

dt^ [dtj ~ ' dr 

Here (r, 0) define the cylindrical coordinates of the particle referred to an origin at the centre of 
mass of the planet. Here ^(r) denotes the gravitational potential due to both the central planet, the 
neighboring satellites and the ring. In addition {Fr,Fo) denote the radial and azimuthal components 
of any additional force F per unit mass respectively. This may arise through internal interactions 
between particles that might lead to an effective pressure and/or viscosity. 

We introduce a Lagrangian description in which the system is supposed to be perturbed from 
an axisymmetric state in which particles are in circular motion with coordinates such that r = vq, 
9 = 9q = r2(ro)t + (3q. Here ro is the fixed radius of the particle concerned, r2(ro) is the angular 
velocity and Pq is a phase factor labeling each particle. In keeping with a Lagrangian description 
(fOiPo) are conserved quantities for a particular particle and so may be used to label it. 

In order to describe the system when it is perturbed from the axisymmetric state we introduce 
the components of the Lagrangian displacement ^ = (Cn^e)- These are such that the coordinates 
of each particle satisfy: 

r = ro + ^r, (3) 
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and 

ro{O-eo)=^0. (4) 

To obtain equations for and we take variations of Eq.'s (1) and (2). We do this by 
applying the Lagrangian difference operator, A, as defined by Lebovitz (1961) to both sides of 
Eq.'s (1) and (2). For a given quantity Q, the variation A(Q) is defined by: 

A(Q) = g(ro + 0-Qo(ro), (5) 

where Q and Qo are the values of the given physical quantity in the perturbed and unperturbed 
flow respectively. In contrast, the Eulerian difference operator is defined as: 

= Q(ro)-Qo(ro). (6) 
Thus, to first order, they are related through: 

A = (5 + ^ . V (7) 
which gives the linear form of the Lagrangian difference operator. 



2.1. Equations for the Lagrangian displacement 

Following Shu et al. (1985) we assume that the components of the displacement are small 
enough that they can be treated as linear in the sense that |^/ro| << 1. On the other hand 
the radial gradient of the radial displacement may be large so that \{d^r/drQ)\ may be of order 
unity. The significance of these assumptions is that although the ring eccentricity is assumed to 

be everywhere small, the ring surface density perturbation induced by it may be of order unity. 
Adopting them enables us to perform the variation in the accelerations using the linear form of 
the difference operator as described above, wherever radial gradients are not involved. These then 
satisfy: 

Here the potential due to the satellite, ■0s j and that due to the self-gravity of the ring, V^SG) are 
included in -0'. Thus 0' = "ipsG + V'S- The quantities = A(i<V), fe = A(i<g) denote the variational 
components of the force per unit mass due to particle interactions. The full non linear Lagrangian 
variation is retained for 0' and F as these may involve the density variation. Contributions coming 
from the variation of the central planet potential are included on the left hand sides of Eq. (8) and 
Eq. (9). 
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2.2. Surface density perturbation and Lagrangian time derivatives 

We suppose the ring particles to be in eccentric orbits and combine to form a globally eccentric 
ring. This is described using a surface density distribution S(r, 9) and eccentricity distribution e(r). 
We also consider there to be an axisymmetric reference state for which e(r) = ^r/^o from which 
we can regard the eccentric ring as being the result of a perturbation. The perturbation of the 
surface density is of the form: 

S(r,^,t) ^ S(r,0) + S'(r,^,i). (10) 

For linear perturbations S' oc cos/sin(m0), where the azimuthal mode number, m = 1. The 
eccentric ring can be thought of as being predominantly in a mode with azimuthal mode number 
m = 1. In practice we may assume |e| << 1. 

We further remark that the convective derivative d /dt is taken following the fluid motion. In 
the approximation scheme used here in which the displacements and hence Lagrangian velocity 
perturbations are small, we may replace the fluid motion by its unperturbed value. Then for any 
quantity Q 

^ = ^ + ^^^. (11) 

dt dt dOo ^ ' 

Similarly 

df^ ~ dt^^ del dOodf ^ > 



3. Forcing of the m=l (eccentric) mode 

We here consider the forcing of the m = 1 mode which causes the ring to become eccentric. The 
forcing is assumed to be due to an external satellite with mass Ms- With an aim of application 
to the Titan ringlet around Saturn we consider that the perturbing potential acting on the ring is 
stationary in a frame rotating with the mean orbital rotation rate when viewed from an inertial 
frame. In terms of the ring dynamics this pattern rotates at a low frequency, 17^, such that VLp «Vt. 
A free m = 1 mode that is most easily excited is one that has a global structure in the ring and has 
a pattern that precesses at a rate comparable to Q,p. When this precession rate is equal to Vtp. there 
is the possibility of resonance and a large response to forcing. The quantity sets the natural 
time scale for variations associated with the m = 1 mode for the problem on hand. 

Accordingly: 

'«"f^) (13) 



dt \d9o. 

Recalling that the left hand side of Eq. (9) approximated by the linearized form, gives for the 
azimuthal component of the displacement: 

^ + 2n^r = Qeo (14) 
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where the quantity Qg^^ is defined by: 



ot o9q r o6q 

Using (13) gives the adequate approximation: 

Wc comment that the motion is dominated by the central mass and to the lowest order Keplerian, 
This means that the m = 1 component of the displacement satisfies (Shu et al. 1985): 



<r- (17) 



Furthermore (14) tells us to lowest order in which Qg^ and ^ may be neglected that: 
which applies to Keplerian orbits with small eccentricity. 

Also using (14) and Eq. (8) one finds that the m = 1 component of the displacement satisfies: 

Here the square of the epicyclic frequency is given by: 

= ^Mi^. (20) 
ro dro 



4. The condition for a steady state response in the rotating frame 

The m = 1 mode responsible for the ring eccentricity has a constant and very small pattern speed 
as viewed in the inertial frame. This means that individual ring particles appear to be in elliptic 
orbits that process at the same rate. In oder to achieve this the internal and external forces acting 
in the mode have to satisfy a constraint that can be view as a non-linear dispersion relation. Our 
treatment again follows that of Shu et al. (1985) who provided such a relationship for density waves 
in Saturn's rings. Except here we consider a density wave comprising a global normal mode rather 
than a forced propagating wave. 

Eq. (19) can also be written in the form: 

^ + &>^' = fr-^ + 2nQe, (21) 
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We now use an angle that is fixed with respect to a coordinate system rotating at the pattern 
angular frequency Op, namely (po = 9q — Qp t. The radial displacement is taken to be of the form 
= -^{ro) cos((/)o)- Following Shu et al. (1985) we note that as the time dependence is contained 
within (f)Q, only depends on ro and (f)Q. 

Multiplying Eq. (21) by cos(^o) and integrating over (pQ, we obtain: 

^2 \ 1 



1 Mro) = TTT-Fr^i^cr + 9D{ro) 



where: 



r>27r 
'0 

10 



1 r'^'' 

+9R{ro) + ^J^ 2 nQgoCos{(Po)d(l)o), (22) 
1 /■^'^ 

^cr = ^j^ /rCOs(0o)#O, (23) 
1 /"^'^ 

9D{ro) = cos(^o) 



2- 

dr 



d(t>o (24) 



and 



1 f'^'^ dips 
c?R(ro) = -— J cos((/>o)-^d(?!)o (25) 

gives the forcing due to the satellite potential which is here considered to be responsible for the 
excitation of the m = 1 mode. 

With the use of equation (16), the last term in Eq. (22) can be re-written and after an inte- 
gration it reads as: 

l-J^\nQg cos((^o)#o = -2 {F^ + gr), (26) 



where: 



and 



l'2n 
10 

10 



^c0 = ^j^ /0sin(,^o)#o, (27) 



Eq. (22) then becomes: 



where 



9T = -^j^ sm<^o-^#o. (28) 
' ^ ll A{r,) = (29) 



and 



2\{n-^pf J ^ {n-^p)^ 

9int = {Far - 2 F^e) + QD (30) 

9s = 9r-'^ 9t- (31) 
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Given that k = ft — uiprec^ where a;prec('"o) is the local radius dependent precession frequency and 
assuming that fip << O and Uprec << Eq. (29) can be approximated to first order in Sip and 

^prec ^* 

(Op - UJprec) A{^o) = (32) 

Eq. (32) provides a condition to be satisfied by the excited m = 1 mode amplitude It balances 
the ring self gravity, internal collisional terms and satellite forcing. For simplicity we shall neglect 
collisional effects below and thus replace gint by go- Note further that for a thin ring of the type 
considered here, may be taken as constant in (32) and evaluated at the ring centre from now on. 



4.1. The self-gravity term 

In order to calculate gn we follow Shu et al. (1985). As radial variations are much more rapid 
than azimuthal ones, the local self-gravity at ro is canonically and adequately approximated to be 
that due to an infinite plane sheet of radial width Ar = r2 — ri , where ri and r2 are the inner and 
outer bounding radii of the unperturbed ring respectively. Thus: 



, . / - - /" 7^^-' (33) 
\ dr J {r-r') 

where G is the gravitational constant, r = ro + and r' = rg + where 

& = Cr(r-o) = A{ro) cos{^o) (34) 

and 



Possible singularities in the integrand are dealt with by evaluating the integral in a principal 
value sense. In the planar limit, we identify the ring eccentricity as e(ro) = 2A{ro)/{ri +r2). Using 
the tight-winding approximation we have: 

T,{r')dr' = S(ro)(iro (35) 

which represents conservation of mass. Since ro/r 1, we have: 

f^U2G r f^^'^\ (36) 
\ dr ) ro + C.-r^-e 

We can re-write Eq. (36) in terms of the eccentricity gradient, q: 



A{ro)-A{r',) 
ro - rA 



(37) 
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Then after integrating over ^, we obtain (see also Shu et al. 1985): 



where: 



Notice that the integrand in Eq. (38) presents a singularity to be handled in a principal value sense. 
This can lead to practical complications near ring edges. 



5. Two integral relations 

The practical problem is to solve equation (32) for the response to the forcing by the external 
satellite. This is equivalent to calculating the forced eccentricity. To do this requires an accurate 
specification of the ring surface density profile which may not be available. Instead one may derive 
two integral relations which contain complete in formation about the response when it is a linear 
function of radius. That is equivalent to assuming the constancy of q defined above which is a 
frequently adopted approximation in planetary ring dynamics (eg. Goldreich and Tremaine 1979, 
Borderies et al. 1983, Shu et al. 1985, Chiang and Goldreich 2000). 

A justification for this is that normally (as here) one considers the case when ring self-gravity is 
strong enough to balance differential precession. In that case the ring precesses at a uniform rate 
similar to a rigid body. When this process is effective, strong self-gravity precludes short wavelength 
displacements so that approximating the induced displacement as a linear function of the distance 
to the ring centre is reasonable as long as the eccentricity response is not too large. 

For example, a strict resonance between the pattern speed and the uniform precession frequency of 
the ring might result in large eccentricities being excited for which the dependence of the precession 
frequency on eccentricity should not be neglected. We note that if > \{Aa/ujprec) x iduJprec/dr)\, 
Aa being the ring semi-major axis width, that effect becomes comparable to that due to differential 
precession and should not be neglected -as we have done- in comparison to that. However, this 
situation is not encountered for the application considered here. 



5.1. The first relation 



To obtain this we take equation (32) multiply by S(ro) and integrate over the ring to obtain 
J S(ro)r2(r2p - ujprec) A{ro)dro - J T,{ro)gDdro = J T,{ro)gsdrQ (40) 
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Prom equation (38) the second integral on the left hand side of equation (40) is zero. The first 
relation thus simplifies to become 

/ S(ro)0(Op 

^prec 

) Airo)dro = J S(ro)55dro (41) 



5.2. The Second relation 

To obtain this we take equation (32) multiply by T,(ro)A{ro) and integrate over the ring to 
obtain 

J ^{roM^p - cuprec) A'^{ro)dro - J ^{ro)A{ro)gDdro = J J:iro)A{ro)gsdro (42) 

In this case we use equation (38) to evaluate the second term of equation (42), which, after 
making use of the symmetry properties of the integral, gives: 

^{ro)A{ro)gDdro = G H H Mj:^r',Mro) (-^(^o) - ^(rj,))^ ^^.^^^ ^^3^ 



I 



which is positive definite. This may also be written entirely in terms of q as 



[ ^{ro)A{ro)gDdro = G H H I{q)q^{r'o)^ro)dr'odro. (44) 



If we now specialize to the case when q is constant we accordingly write ^(ro) = Ac + qx ss & 
linear function of radius. Here Ac is constant and x = — rc measures the radial coordinate 

relative to the centre of mass of the unperturbed ring assumed slender. The eccentricity at the ring 
centre satisfies is e = \Ac\/rc, while the sign of Ac determines whether pcriccntre is in the direction 
00 = 0, or TT if negative or positive respectively. This means that by definition 

j x^dx = 0, (45) 

Jri 

where the domain of integration for the above and similar integrals below is the extent of the ring. 
We further adopt a linear form for the precession frequency, thus 

t^prec = ^precfi + '^'precfi^^ (4^) 

where ct;prec,o and oj'p^cc,o represent the precession frequency and its derivative evaluated at x = 0. 
Consistent with the approximations made here we may also assume the forcing potential term gg 
is constant throughout the ring. 

The first relation (41) then gives a first relation between Ac and q in the form 

j T,Q{Qp - LOprecfi) Acdx - J 'EQuj'pj.ccfiqx'^dx = J T,gs.dx (47) 
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The second relation similarly leads to a second which takes the form 

J T,Q{^p — iOprec,o) (^c + <f'X^)dx ~ j '^^'^precfi^cQX^ dx = 

Y.^J^^^^,^(^x^dx + G / I{q)qT,{x')'E{x)dx'dx + / ^AcQsdx. (48) 

J Xl J Xl J 

We now have two relations which enable the forced response to be calculated under the assumption 
of constant q. These are A,., whose magnitude, when divided by the radius at the ring centre gives 
the eccentricity at the ring centre and q itself which is the product of the central radius and the 
eccentricity gradient. Note that implicit in the thin ring approximation is the requirement that the 
magnitude of the eccentricity gradient significantly exceeds the ratio of the eccentricity to radius. 



6. The relation between q and central eccentricity 

If we multiply the first relation (47) by Ac and subtract it from the second relation (48), the 
terms involving satellite forcing cancel out and we get a relation between Ac and q which (recalling 
that q is constant and fl is evaluated in the center of the ring) takes the form. 

q^iClp — UJprec,o) J ^X^dx — ^LOpj.gc^QAc J TiX^dx = 
qnu'p^^cO [^x^dx + G I ' / ' I{q)J:{x')^{x)dx'dx. (49) 
This may be cast in the very simple form 

. ^ g(Op - (jJprec,o) _ q J T^x^dx _ GI{q) (J T.dxf 



6.1. Determination of q 

We may now use the relation between q and Ac specified above, to eliminate Ac in the first 
relation (47) and so obtain q in terms of the satellite forcing. This gives 

/ {Q.p - UJprecof {^P - i^preco) J X^T^dx J T^x'^dx \ 



Hrec,o)^ ^'precfi J x'^T.dx J T,dx 



(J^P - ujprec,o)GIiq) (/ T,dx)^ gs 



^Hreco)^ / ^x'^dx ^i^precfi ' 



(51) 



Thus equations (50) and (51) give the response parameters Ac and q directly in terms of the 
external forcing. 
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6.2. The linear regime 

When the response is in the linear regime q is small and I{q) = q/2. In this case the right hand 
side of (51) is proportional to q and we have 



where 



_ (r^P - U^preco)^ _ J Sx^dx _ (Op - U;prec,o) J X^T.dx 



{Up - UJprec,o)G (/ Sdx)^ 



The response is singular when ai = 0. Regarding this as an equation for Op — ujprec,Oj we have 
a quadratic with two real roots indicating a singular response for certain pattern speeds. However, 
the relation (51) is in fact a nonlinear one (through the functional form of I{q)) to determine q and 
the nonlinearity present can remove such singularities. This becomes apparent when we expand 
to the next highest order in q. Then one obtains a cubic equation for q that can always be solved 
because the coefficients of q^ and q never vanish simultaneously. This cubic takes the form 

_ 3 3(Op - U}prec,o)G (/ Sdx)^ _ gs 

However, as indicated above we must be cautious about using the above determinations of Ac and 
q when the eccentricity response is large because then the assumption of constant q and the neglect 
of the dependence of externally induced orbital precession on the eccentricity may no longer be 
valid. 



7. Application to the Titan ringlet 

We here consider the Titan ringlet for which the ring precession rate is close to the orbital 
frequency of the satellite Titan. This ringlet is therefore a candidate for having its eccentricity 
forced by Titan. We now consider the forcing potential. 



7.1. The satellite potential 

For the low frequency m = 1 forcing considered here the perturbing effect of the satellite arises 
through 

'' = -^1 (^°^('^°)^-2«^-('^o)^||)#o (55) 
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For the satellite we neglect the orbital eccentricity and expand tpg to leading order in r/05, 05 
being the semi-major axis of the satellite orbit. For the forcing considered here, recalling that in 
the frame rotating with the orbital frequency 0,p, only secular terms oc cos(^o) are significant, we 
obtain, including the indirect potential, to leading order 

/ ^GMsr^ . 

■08 = — cos((?!)o). (56) 

Then evaluating at the ring center r = Tc = a, we obtain 

We are now ready to apply equations (50) and (51) to the Titan ringlet. 



7.2. Parameterizing the dynamical model 

We define the dimensionless parameter rj through 

rip - Uprec,0 = r] (^'precO (58) 

This defines the resonance where the precession rate induced by the planet and Titan match in the 
approximation that the former can be represented by a first order Taylor expansion about the ring 
centre. It gives the resonance location at a distance rjAa from the ring centre. 

Wc also find it convenient to define the dimensionless quantities 7„ related to the ring surface 
density profile through 

The location of the center of mass of the ring,rc, is the origin of the coordinate x (equation (45)). 
Then, inside the ring we have |x| < Aa -the extreme case is when the mass is concentrated at a 
ring edge. Thus, it is verified that |7„| < 1. We also notice that 7n > when n is even and, when 
n is odd, 7^ = if the surface density is symmetric. 

Further wc introduce the parameter $ which measures the importance of self-gravity with 
respect to differential precession: 

where the mean surface density is defined as: 

notice that since the precession frequency decreases with distance (i.e. w^^g^ g < then $ is 
defined as a negative quantity. 



-14- 



In terms of these dimensionless parameterizations, equations (50) and (51) can be reduced to 

Ac = qAa 02, (62) 

where ^ 

a2 = ri- — , (63) 

72 72 

and 

7.3. Estimation of the surface density and the location of the secular resonance 

The presently available observations of the Titan ringlet do not enalble one to determine all 
its dynamical and physical parameters. Optical depth profiles give the most accurate description 
of the ring, thus the surface density can only be inferred by making assumptions on the physical 
properties of the ring particles or through dynamical models which, like the present one, rely 
on various assumptions (see for example Goldreich and Tremaine 1979). Moreover, there are 
considerable uncertainties in the values of the multipole moments of Saturn (for up-to-date values 
see the JPL-Solar System Dynamics website: http://ssd.jpl.nasa.gov/sat-gravity. html). In fact, the 
uncertainties in the precession frequency of the ring due to the uncertainties in multipole moments 
of Saturn, imply an error in the location of the secular resonance that is of the order of the width 
of the ring. 

Using our model we shall attempt to put some constraints on the value of $ and so, on the 
mean surface density , E, as well as on the resonance location parameter, r), for the Titan ringlet. 

We shall adopt the following values: a = 77871A;m, Aa = 25km, e = 2.6x10"^, 6e = (1.4±.4) x 10~^, 
(Porco et al 1984). Thus q = ade/Aa « 0.44 ± 0.18. 

The closest approach between the satellite and the ring occurs at apoapsc (Porco ct al. 1984). Then, 
according to the definitions of the satellite potential (equation (56)) and of the radial displacement 
(equation (34)) adopted here, when ^ = 0, > 0. Thus, for the Titan ringlet we have Aj.{ro) > 
as well as > 0. 

Now, we shall rewrite equations (62) and (64) as: 

as = ?? 02 - 72 (65) 
-$ = (02 -rj) j2 + 7s, (66) 



where: 



Notice that 03 < since uj' c < ^ 55 > 0- 
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We can eliminate r] from equations (65) and (66), to obtain a quadratic expression of $ as a 
function of the form factor 72: 

-$ = _ J_ ^2 ^ (02 - ^) 72 + 73- (68) 
a2 a2 

Similarly, one can eliminate 72 and express $ as a function of rj: 

-$ = -02 77^ + {al + 03) 77 — 0302 + 73- (69) 

One must recall that 77 and 72 are related by equation (65), which acts as a constraint for 77 since 
it must be verified that 72 < 1. 

Note also that 02 can be expressed as a relative eccentricity, 

through the use of equation (62). The quadratic terms in equations (68) and (69) depend on 02- 
Then, if the relative eccentricity is large, is: — $ « 0272 + 73- Whereas if 02 is small, — $ « 
03 77 - 0302 + 73. 

One may expect to find observed values of 72 between those of two critical cases. One extreme 
case for the shape of the surface-density distribution is when it is constant throughout the ring, 

E«/E«(x), 

the other when the whole mass is concentrated at both edges, as: 

(X 5{x - Aa/2) + 5{x + Aa/2), 

were 5 is the Dirac-(5 function. In these cases the values of 72 are given by 72^'' = 1/12 and '^^^ = 1/4 
respectively. 

We can estimate the values of 02 and 03 for the Titan ringlet as: 



02 



(^) « 1.86 (70) 
\5eJ Obs 

03 -1.38, (71) 

where these estimates are obtained from the observed parameters previously quoted and, in equa- 
tion (71), the precession frequency gradient, i^p^ecO' computed using the J2, J a and Jg cof- 
ficicnts associated with the multipole moments of Saturn given in the JPL-SSD website: J2 = 
(16292 ±7) X 10-^ J4 = (931 ±32) x 10~^ Jg = (91 ±32) x 10"^ (Jacobson 2004). We computed 
the precession rate of the ring due to the oblateness of Saturn and the secular perturbation by 
Titan ^ and we obtain a value of a^es = (77846. ± 11.) km, for the location of the 1 : secular 
resonance, where the errors only consider the uncertainties in the values of Saturn's J-coefficients. 



^We use an expansion up to {Rsaf/af' in the oblateness term. For the secular term, we used an approximation 
up to first order in e for the perturbing function (see for example Murray & Dermott 1999), where we have used a 
value of aTitan = 1221900. km for the semimajor axis of the satellite and of rriTi tan — 1345.5 X W'^^gr, for its mass. 
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For the symmetric case, we have plotted the relationships set by equations (65) , (68) and (69) 
for the values of 02 and 03 corresponding to the Titan ringlet. 

We may recall that, if the surface density is symmetric with respect to the center of mass, then 
73 = and that inside the ring |x| < Aa/2 and \r]\ < 0.5. 

Figure (1) shows S vs. 72 -where E is related to $ by equation (60). It can be seen that 
the scale of the mean surface density is tens of grammes per cubic centimetre, in particular when 
72^^ < 72 < 72^^- symmetric case, the agreement with the solution quoted in Porco et 

al. (1984) occurs at a value of 72 ~ 72^'*, however, in that case the numerical values of various 
parameters are different. 

In figure (2) we plot t] vs. 72. This figure gives, for a symmetric profile, the location of the 
secular resonance as a function of the form of the surface density, such that the ring precesses 
rigidly when the effects of inter-particle collisions are negligible. The allowed values of r), set by 
equation (65) and 72 < 1, imply that, given the condition of solid precession, the location of the 
secular resonance lies outside the ring (I77I > 0.5) for the symmetric case with 72^'* < 72 < 72^^- 
We may recall that if r/ ~ 0, then the assumptions made to obtain these relationships break down, 
because the eccentricities resonantly-excited may become very large. Thus, the rather large values 
of T] obtained assure us that the first order approximation adopted here is consistent with the case. 
However, as mentioned earlier, the uncertainties involved are large. 

Finally, in figure (3) we plot S as a function of tj. This figure presents the required values of the 
mean surface-density to produce the balance as a function of the location of the secular resonance. 
Naturally, the result reflects the particular choice of parameters adopted. All the solutions obtained 
require rj < 0, hence the the center of mass of the ring must be located interior (i.e. towards the 
satellite) with respect to the secular resonance (see equation (58)). 

If the surface density profile were not to be symmetric then these estimates can be somewhat 
altered. By definition, the magnitude of 73 is bounded by that of 72 as: 

I73I < 72 < 72- 

In figure (1) we also show S as a function of 72, for the extreme cases in which I73I =72. 

On the other hand, we do not expect the effect of the collisional terms to be negligible, as 
assumed in these calculations. We will explore their exact significance in a forthcoming paper. We 
will discuss in the next section how physical particle interactions could affect the results obtained 
here. 



8. Discussion 



In this paper we have used a model of a thin slender self-gravitating ring in orbit about a dom- 
inant central mass to produce relationships between its physical parameters. These relationships 
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can then be applied to real systems to make predictions about their physical and/or dynamical 
properties. 

We view the ring in uniform precession as sustaining a global non-axisymmetric m = 1 mode 
of oscillation. We have considered particularly the case in which the precession frequency is in a 
secular resonance with the orbital frequency of an external satellite, which is the case of the Titan 
ringlet of Saturn. 

A condition for the ring to be able to maintain a m = 1 mode with single slow pattern 
speed is obtained. It can be expressed as an integral condition for the pattern or normal mode to 
precess at a uniform rate (Eq. (29)) that requires the correct balance between differential precession, 
self-gravity, secular-resonant forcing and collisional effects. Prom this condition we have obtained 
two relationships (equations (41) and (42)) which combine different observable parameters of the 
ringlets. 

We argued that the scale of the perturbations is of the order of the width of the ring, because 

we are interested in cases in which the differential precession is compensated mainly by self gravity 
and the eccentricity is small. Hence, a first order expansion is a good approximation and so, we 
take the eccentricity gradient as constant, which, indeed, is the case in real systems (Graps et al. 
1995). 

In the linear case, the two expressions obtained give simple relationships (equations (47) 
and (48)) between the eccentricity at the center, ec, the eccentricity gradient, q, the location 
of the secular resonance, rj and the mean value of the surface density S, and its form, expressed by 
the 7n's. The relationships can be further simplified by assuming that the system is also linear in 
q (equations (50) and (52)). 

Finally we have applied this model to the Titan ringlet. The scale of the surface density ob- 
tained is consistent with previous estimations (Porco et al. 1985), set at the order of ~ 10 gr cm~'^. 
Our symmetric solutions imply that the secular resonance is outside the ring (see figure (3)). How- 
ever, the distance between the ring and the location of the secular resonance obtained is smaller 
than the total uncertainty in the location of the later, given the errors in the determination of the 
multipole moments of Saturn. 

Moreover, additional physics may need to be considered. The fact that all the observed systems 
have shown a positive value of q is an indication that the differential precession is mainly balanced 
by self-gravity (Borderies et al 1983, Papaloizou and Melita 2004). However, it is also very likely 
that physical collisions play an important role in setting up the rigid precession, particularly at the 
edges, where the relative motion is larger (see for example Chiang and Goldreich 2001). It is also 
believed that the real systems are close-packed at their narrowest point -as it seems to be indicated 
by the observations of the Uranian rings (see for example French et al. 1984). The narrowest 
points of the observed eccentric ringlets occur at periapse, since q > 0. Thus, the collisional 
contribution can be estimated to be impulsive, i.e. entirely concentrated at the pinch. At periapse, 
the geometry of the collisions at the edges is such that the coUisionally produced orbital phase-shift 
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produced reinforces the differential-precession induced by the oblateness of the central planet (see 
for example Moulton 1935). Thus, if a considerable collisional contribution arises mainly at the 
pinch, self-gravity must counteract a greater differential precession. If that is the case, the values 
of the masses of the real eccentric ringlets will be larger than the estimates produced here. In 
fact, the masses of eccentric ringlets estimated from observations have always turned out to be 
greater than the theoretical estimates obtained with models that consider only the self-gravity of 
the ringlet (sec for example Mosqueira & Estrada 2001). It is straightforward to include collisional 
effects when they are assumed to act impulsively at pericentre. -more terms are to be retained 
when considering the two conditions for rigid precession given by (equations (40) and (42)). We 
can gain some knowledge of these collisional effects if an observed profile is available -from which 
we can extract density estimates. We will present a more detailed analysis of some known eccentric 
ringlets in a forthcoming article. 
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Fig. 1. — The estimated mean surface density, S, as a function of 72 for a symmetric profile (73 = 0). 
Notice the value of S when 72 is between the critical values 72"'^'' = 1/12 and 72^'' = 1/4 (indicated 
by the markers). The dotted lines correspond to the non-symmetric case where I73I = 72. 
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Fig. 2. — The location of the resonance r], as a function of 72 (73 = 0) that gives the correct 
balance enabling the rigid precession of the ring, when the collisional terms are neglected. Notice 
that solutions with 7^^^ < 7 < 7^^\ correspond to resonance locations outside the ring (I77I > 0.5). 
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Fig. 3. — The estimated mean surface density, S, as a function of r/ (73 = 0). A negative value of 
rj implies that Op > u)prec,o (equation (58)), i.e. the center of mass of the ring is exterior -from the 
planet- to the location of secular resonance. 
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